Measuring Distance Ratios with CMB-Galaxy Lensing Cross-correlations 
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We propose a method for cosmographic measurements by combining gravitational lensing of the 
cosmic microwave background (CMB) with cosmic shear surveys. We cross-correlate the galaxy 
counts in the lens plane with two different source planes: the CMB at z ~ 1100 and galaxies at an 
intermediate redshift. The ratio of the galaxy count/CMB lensing cross-correlation to the galaxy 
count /galaxy lensing cross correlation is shown to be a purely geometric quantity, depending only on 
the distribution function of the source galaxies. By combining Planck, ADEPT and LSST the ratio 
can be measured to ~ 4% accuracy, whereas a future polarization based experiment Uke CMBPOL 
can make a more precise (~ 1%) measurement. For cosmological models where the curvature and 
the equation of state parameter are allowed to vary, the direction of degeneracy defined by the 
measurement of this ratio is different from that traced out by Baryon Acoustic Oscillation (BAO) 
measurements. Combining this method with the stacked cluster mass reconstruction cosmography 
technique as proposed by Hu, Holz and Vale (2007), the uncertainty in the ratio can be further 
reduced, improving the constraints on cosmological parameters. We also study the implications 
of the lensing-ratio measurement for early dark energy models, in context of the parametrization 
proposed by Doran and Robbers (2006). For models which are degenerate with respect to the CMB, 
we find both BAO and lensing-ratio measurements to be insensitive to the early component of the 
dark energy density. 



I. INTRODUCTION 

Weak gravitational lensing of the cosmic microwave 
background (CMB) (see [1] for a review) provides us with 
a unique opportunity to study the large scale distribu- 
tion of dark matter in the universe out to much greater 
distances than accessible through conventional galaxy- 
lensing studies. It has been shown [2] that by studying 
the gravitational lensing of galaxies in different redshift 
slices by the same foreground structures, the geometry of 
the universe and eventually, dark energy evolution may 
be constrained — a method known as cross-correlation 
cosmography. In this paper, we propose a similar method 
in which we treat the CMB as one of the background 
slices. This not only provides an extremely well stan- 
dardized distance to compare other distances to, but also 
incorporates the longest possible distance in the ratio of 
distances probed by this method, making it a more sensi- 
tive probe of cosmological parameters than ratios involv- 
ing distances restricted to galaxy surveys. Recently, Hu 
et al. [3] have proposed a method for measuring the same 
lensing-ratio by comparing the convergence profile of the 
a cluster reconstructed via background galaxy shear with 
that reconstructed via CMB lensing, and then stacking 
several clusters to improve the precision of the measure- 
ment. The method we propose here depend on cross- 
correlations rather than reconstruction of convergence of 
individual objects, and will have different systematics. 
As such, this method is a powerful complement to the 
cluster-lensing based method and the two may be com- 



bined to obtain more precise measurements of the ratio. 



II. LENSING RATIO: THE KEY OBSERVABLE 

Cosmological weak lensing effects are conveniently en- 
coded in the effective convergence field, which is defined 
as a weighted projection of the matter overdensities 5 [4], 
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where dA{r]) is the comoving angular diameter distance 
corresponding to the comoving distance rj. Here, 0(77) is 
the scale factor, while fim and Hq represent the present 
values of the matter density parameter and the Hubble 
parameter, respectively. The quantity 17(77) represents the 
fact that sources are distributed in comoving distance 
with a normalized distribution function Wb- Since the 
CMB photons all come from nearly the same cosmologi- 
cal distance, we can approximate the source distribution 
function as, Wb{ri) ~ S{r] — 770), giving 
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where 770 is the comoving distance to the last scatter- 
ing surface. We will denote the same quantity for a 
background galaxy population with redshift distribution 
Pg{z)dz = Wb{ri)dri, with the symbol ggediv)- 
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We also consider a suitable foreground population as a 
tracer of large-scale structure. The projected fractional 
overdensity of the tracers can be written as, 



S(n) = J diiWf{r])6g{7jh, 
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where 5g represents the fractional tracer overdensity and 
Wf is the normalized tracer distribution function in co- 
moving distance. We assume that the Fourier modes of 

the tracer overdensity field arc related to those of the 
underlying matter density field via a scale and redshift 
dependent bias factor, so that 5g(k, 77) = b{k,r])S{'k,r]). 
If we cross-correlate the tracer overdensity map with the 
convergence field, we obtain the cross power spectrum, 

C^^ = \n^Hlj drjbMWfiv)^P{£,ri), (5) 

where we have used the Limber approximation and the 
orthogonality of spherical harmonics. We have also in- 
troduced the shorthand notation, be{ri) = b{-^,ri). 

Now, we will introduce two separate cross-correlation 
measures involving the foreground tracer population. 
First, we consider the case for the CMB as the back- 
ground source. By constructing estimators out of 
quadratic combinations of CMB fields (temperature and 
polarization), it possible to obtain a noisy reconstruction 
of the convergence field (note that the actual observable 
in this case is the deflection field) out to the last scatter- 
ing surface [5, 6], which we denote as Kcmb- The power 
spectrum of the noise in the reconstruction, N^cmbkcmb ^ 
can be estimated knowing the specifications for the CMB 
experiment. The cross-correlation of the reconstructed 
convergence field from the lensed CMB with the fore- 
ground tracer, gives the signal, 

^ (6) 

where we have used the source distribution kernel fjcMu 
appropriate for the CMB being the background source. 

Next, we consider the case for the weak lensing of back- 
ground galaxies. The relevant observable in this case is 
the traceless symmetric shear field on the sky, the mea- 
surement of which allows a noisy reconstruction of the 
convergence field appropriate to the background galaxy 
distribution, Kg^i. In this case, the noise is primarily due 
to intrinsic ellipticity of the background galaxies and has 
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the spectrum, jv^''^''"'^''' ~ (Tint) /" where {■jf^f) 
and n is the number of background galaxies per steradian 
[7]. If we cross correlate this convergence field with the 
foreground tracers, we find the signal, 

C;-^ = InmH'oJ dvbeiv)Wf{vf-^P{^,v), 

(7) 

where we have used the source distribution kernel, g'gai 
appropriate for background galaxies. 



If the foreground distribution is narrow in redshift so 
that it can be approximated by a delta function, Wf (jj) ~ 
5{r]—r]f), then the ratio of the above two cross-correlation 
measures, which we call the lensing-ratio, reduces to. 
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which is simply the geometrical ratio of the source dis- 
tribution kernels. If the background galaxy distribution, 
too, is sufficiently narrow in redshift around z = z^al, 
this becomes, 
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Note that this is independent of the angular scale, tracer 
bias and the power spectrum. Therefore, measurements 
at several multipoles can be combined to constrain the 
lensing-ratio. Since the distance ratios depend on the 
cosmology and specifically on the dark energy model, this 
can be used to constrain dark energy parameters. 



III. UPCOMING SURVEYS AND A NEW 
PROBE OF DARK ENERGY AND CURVATURE 

Large scale structure surveys, together with precision 
measurements of the CMB anisotropics have already pro- 
vided us with a wealth of knowledge about the geometry, 
evolution and composition of the Universe. In the com- 
ing decade, Cosmologists will carry out even larger scale 
galaxy and lensing surveys and produce higher resolution 
CMB maps. We consider a combination of three exper- 
iments in order to assess how well the lensing-ratio can 
be measured in such future surveys. We consider the 
redshift slice of foreground tracers (lenses) to be drawn 
from an ADEPT-like [17] large scale structure survey and 
the background (source) galaxies taken from an LSST- 
like [18] weak lensing experiment. For the CMB lensing 
measurements, we consider the upcoming Planck mission 
as well as a prospective polarization-based mission like 
CMBPOL. 

The foreground galaxy slice is taken as a step func- 
tion in the redshift range (0.8,0.9) with 350 galaxies per 
square degree. The source galaxies are also assumed to 
be distributed imiformly in redshift, between z = 1.2 and 
1.6 with a number density of 40 galaxies per square ar- 
cmin. We model Planck to be a 7' FWHM instrument 
with temperature and polarization sensitivities of 28 and 
57 /zK-arcmin, respectively. For CMBPOL, we adopt a 
3' beam FWHM and temperature and polarization sensi- 
tivities of 1 and 1.4 ^K-arcmin, respectively. Wc assume 
that both CMB experiments cover 65% of the sky and 
all cross-correlations are performed over the same area. 
For calculations performed here we assumed a WMAP 5- 
year normalized ACDM cosmology with f2b/i^ = 0.0227, 
ncDuh"^ = 0.1099, Qa = 0.742, r = 0.087, = 0.963 
and As = 2.41 x 10"^. 
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FIG. 1: Cross power spectra, the ratio of which is being stud- 
ied (cf. equation 8). Also shown are predicted la errors in 
uniform bins of size A£ = 30. For the "^"^ case, the 
outer (lighter) errors correspond to lensing reconstruction us- 
ing temperature and polarization with Planck, while the inner 
(darker) ones correspond to the same for CMBPOL. 



In Fig. 1, we display the two cross power spectra ap- 
pearing in the defining equation (8) of the lensing ratio 
along with binned uncertainties predicted from the ex- 
perimental specifications. 

The error on the ratio can be obtained as follows. We 
begin by defining the log-likelihood, 
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where, Zf — 
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We compute the vari- 



ance of Zi at the value tq of r computed in the fiducial 
cosmology, 

1 



a'iZi 



i2e + 



c, 



J.2 ( ^'^gal^gal^ES 



kcmbS\2 



^/^CMB/^gal /^SS 1 /^KCMrS /-^/^galS 



(11) 



where. 
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include the noise power spectra. The Poisson noise for 
the foreground tracer is taken as Nf^ = 1/n/. Then 
maximum likelihood estimate for the ratio is then ob- 
tained by solving d'x^{r)/dr = to be. 
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Now, we can estimate the error on r as, 
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Various auto, cross and noise power spectra that enter 
the calculation of the error on r are shown in Fig. 2. 
The above figures borne out the expected feature that 
the noise power spectrum in the lensing reconstruction 
is the largest source of uncertainty that propagates into 
the error on r. 
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FIG. 2: Various power spectra that enter the calculation of 
the error on the lensing ratio (cf. eq. 13). Each of the noise 
power spectra has been plotted with the same line style as its 
corresponding signal power spectrum and labeled as A'^ . The 
noise spectrum for the CMB lensing reconstruction has been 
indicated both for Planck and CMBPOL. 



The estimated errors on r are shown in Table I. For 
Planck, we find that the lensing-ratio can be estimated 
to ~ 4% while with CMBPOL a ^ 1% measurement is 
possible. 



IV. PARAMETER CONSTRAINTS 

For Planck priors, improvements on cosmological pa- 
rameter constraints upon adding the lensing-ratio to the 
primary CMB observables become appreciable when the 
error on the ratio decreases below 10% [3]. It is inter- 
esting to note here that the method for estimating the 
lensing-ratio proposed by Hu, Holz and Vale (2007) [3], 
which relies on cluster mass reconstruction can be further 
improved with the maximum likelihood based estimator 
proposed by Yoo and Zaldarriaga [8] and can comple- 
ment the method proposed here. By combining the two 
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FIG. 3: Left Panel: Improvements of constrains in the Qk — f^A plane for a vacuum energy model with Planck by adding a 1% 
measurement of the lensing-ratio. The outer solid contour is the 68% confidence interval from primary CMB alone while the 
inner solid contour is the same after adding the lensing-ratio. The dotted contours have the same interpretation but represent 
the case where information from lensing extraction has been added to the CMB Fisher matrix. Right Panel: Same as left, but 
for the w — JIa plane, assuming flatness. 



Experiment Type (S/N)"°== Ar/r(%) 



Planck 


POL 


25.8 


3.8 




TT 


23.3 


4.2 


CMBPOL 


POL 


102.6 


1.0 




TT 


84.5 


1.2 



TABLE I: Predictions for the cross-correlation studies de- 
scribed in the text with foreground galaxies from ADEPT, 
background galaxies from LSST and different CMB experi- 
ments. The quantity (S/N)'^'^"'''' represents the signal-to-noise 
ratio in the estimation of the cross correlation between the 
foreground tracer density with CMB lensing. The last col- 
umn shows percentage error in the lensing-ratio estimator, r 
of equation (8). We show the prediction for both temperature 
based (TT) and polarization based (POL) reconstruction of 
the deflection field from the lensed CMB. 



methods for the same redshift slices, it may be possible 
to reduce the uncertainty in the lensing-ratio to percent 
or sub-percent levels. 

In order to assess how a percent-level measurement of 
the ratio will help constrain a set of cosmological param- 
eters {pi} in conjunction with the CMB experiments, we 
define a Fisher matrix for the lensing-ratio. 



dlnr 1 dlnr 
■ ^ _ 

dpi (In r) dpj 



(14) 



and add it to the Fisher matrix from a CMB experiment. 
The error in a parameter is then estimated from the in- 
verse of the combined Fisher matrix as cr(pi) = -y/fF^^. 
We consider two variants of the CMB Fisher Matrix, one 



with only the primordial power spectra and the other 
with the power spectra involving the weak lensing de- 
flection field extracted from CMB lensing measurements 
[9, 10]. We do not consider any foreground contamina- 
tion in any of these. Fig. 3 shows the constraints pre- 
dicted with Planck specifications and a 1% error on the 
lensing-ratio, for minimal extensions to the standard 6- 
parameter model. These constrains are marginalized over 
all other parameters. The constraints on curvature as- 
suming w = —1 and on w assuming flatness, both im- 
prove over the primary CMB case after adding in the ra- 
tio. For CMB with lensing extraction the improvement 
on w is still substantial while that on is marginal. 

Following [9], we also consider a more general 11- 
parameter model with w, fraction of dark matter in 
massive neutrinos f^, the effective number of neutrino 
species, the running of the spectral index and the pri- 
mordial Helium fraction. We adopt a flducial value of 0.1 
eV for the total neutrino mass. The constraints on the 
interesting subspace of parameters are shown in Fig. 4. 
We find in this case that adding in the ratio significantly 
improves the constraints on w and over the primary 
CMB case. In fact, for w and J^a the improvements sur- 
pass those from the lensed CMB Fisher matrix. This 
is particularly interesting because the lensed-CMB-only 
constraints require an estimate of the convergence field 
from the four point function in the lensed CMB itself and 
is more prone to systematics than the cross-correlations 
that enter the ratio calculation. From Fig. 4, it is ap- 
parent that for the CMB Fisher matrix with lensing ex- 
traction the constraint on the neutrino mass is rather 
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FIG. 4: Improvements in the constraints on the interesting subset of parameters in the eleven parameter model involving 
massive neutrinos and free dark-energy equation of state (see text). The interpretations of the contours are same as in Fig. 3 



tight, SO that no further gain is obtained by adding in 
the lensing-ratio. It is important to keep in mind that 
Fisher matrix methods tend to overestimate the error on 
neutrino mass due to non-Gaussianity in the associated 
hkelihood [10]. Therefore, the errors estimated here are 
somewhat higher than those predicted from a full Monte- 
Carlo forecast. 

We next consider models with curvature and a free 
dark energy parameter w. CMB measurements allow for 
a large degenerate valley in the w — ^Ik plane making 
Fisher Matrix results deceptive. We treat this model by 
imposing a strong CMB prior in which we explore the 
parameter space while keeping the high-redshift variables 
ilhh"^, ^mh^ and 9a = rs{zo)/dA{'rio) fixed. For each de- 
generate parameter set, we calculate the value of r as 



well as the spherically averaged baryon acoustic oscilla- 
tion (BAG) distance ratio [11, 12], rs{zd)/ Dy{z), where 
rs is the comoving sound horizon scale at the drag epoch, 
Zd and is an effective distance measure to redshift z. 
We compute it for z = 1.5, a typical median redshift for 
an ADEPT-like survey. The constraints in the w — fife 
plane from each of these methods are shown in Fig 5 for 
a 0.4% measurement of the BAG ratio with ADEPT [13] 
combined with a 1% measurement of the lensing-ratio. 
We find that the degeneracy direction for the BAG is 
quite different from that for the lensing-ratio. Together 
they can put a ~ 0.01% limit on fife and a simultaneous 
10% limit on w, without adding in any other cosmo- 
logical prior. 

Next, we turn to scenarios with early dark energy. In 
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FIG. 5: Constraints (68%) in the w-Q.k plane from BAO and 
lensing-ratio measurements. The dashed hne indicate con- 
straints from the lensing-ratio while the dotted line represents 
the same for the BAO ratio. The solid contour shows the joint 
constraint. Regions outside the contours labeled "excluded" 
are not allowed due to the strong CMB prior (see text). 



again impose the strong CMB prior, and compute the 
obscrvables for each point in the (wo,^^^) space degener- 
ate with respect to the CMB. Since early dark energy 
shifts the comoving sound horizon according to equa- 
tion (17), and we are fixing the values Q.mh? and Q.bh? 
(which fixes the rcdshift of last scatter), equation (18) 
immediately implies that the only way to keep the angu- 
lar scale 9a = rs{zo)/dA{i]o) constant is by varying Vl^ 
such that (Ia scales like ^1 — fi^. Thus, quite counter- 
intuitively, we find that the early value of the dark energy 
indirectly affects low redshift evolution. An unfortunate 
consequence of this is that both BAO and lensing ra- 
tio measurements are rendered insensitive to the value of 
Cl^. As shown by [15], it is possible to have a set of cos- 
mological models degenerate with respect to the CMB 
without having flmh^ strictly constant. For example, 
holding f2ni, and therefore fl'^ constant, it is still possible 
to have nearly indistinguishable CMB power spectra. In 
this case, to preserve 9a, one has to change h accordingly. 
This, again, leads to the shift in low rcdshift distances 
and makes our observables insensitive to Q^. These is- 
sues are discussed in some detail in the Appendix. 



particular, we choose the dark energy parametrization 
proposed by [14], namely, 

where is the present value of the dark energy density 
function, fid (a), and is its asymptotic value at high 
redshift. In this parametrization, the dark energy equa- 
tion of state w{a) has the value wq at present, crosses 
over to w ~ during matter domination and goes to 
w ~ 1/3 in the radiation dominated era. Two relations 
of interest in this model that follow from the definition 
of the Hubble parameter. 



2 (Oma ^ + ^rtt '*) 



1 - naia) 



(16) 



are the scaling, with f2^, of the comoving sound horizon 
at last scattering or the drag epoch, namely. 



(17) 



and the behavior of the comoving angular diameter dis- 
tance which, for flat cosmologies, given by, 
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dz 



[fl^{l + z)^ + Qr{l + zy]y^' 



(18) 



We consider the parameter space spanned by (ti;o,fi^) 
and study how well it can be constrained given measure- 
ments of the lensing ratio, the CMB and the BAO. We 



CONCLUSIONS 



We have proposed a way of measuring a ratio of comov- 
ing angular diameter distances that appear in the lens- 
ing kernels for CMB and galaxy lensing. By combining 
Planck, ADEPT and LSST, it is possible make a percent 
level measurement of this ratio. A polarization based ex- 
periment like CMBPOL has the potential of making a 
more precise measurement. The precision in the mea- 
surement can be potentially increased by combining it 
with cluster mass reconstruction based measurement of 
the same quantity. The ratio is sensitive to late-time ge- 
ometry and composition of the Universe and a percent 
level measurement combined with Planck data can pro- 
vide interesting constraints and consistency checks, inde- 
pendently of other cosmological probes. By choosing the 
CMB as one of the lens planes, this method allows higher 
redshifts to be probed than galaxy-lensing cosmography. 

While the distance ratio is sensitive to late time dark 
energy, we find it to be rather insensitive to early dark 
energy, particularly for the parametrization proposed in 
[14]. As discussed in the Appendix, when a strong CMB 
prior is imposed, the values of low redshift parameters 
shift in conjunction with the asymptotic high redshift 
value of early dark energy in such a manner as to ren- 
der both the BAO ratio and the lensing ratio rather in- 
sensitive probes. This behavior is most likely a specific 
feature of the said parametrization. The effectiveness of 
the lensing ratio as a cosmological tool for a wider class 
of quintessence models remains to be studied. 



7 



APPENDIX: BAO AND LENSING RATIOS AS 
PROBES OF EARLY DARK ENERGY 

Let us choose the fiducial dark energy model as one 
with n'^j = 0.742, Wo = -1 and fl^^ 0.03 with the other 
cosmological parameters set at the values described in 
the body of the paper. Now, if we had wrongly assumed 
^l^ — i.e. a A model, then by fitting the CMB, we 
would find a set of parameters that would keep the angle 
= '''s{zo)/dA{zo) constant. As oc (1 — fl^y^^, we 
would find a model that overestimates the sound horizon 
and hence the comoving angular diameter distance to the 
last scattering surface by the factor (1—0.03)"^/^ ~ 1.015 
i.e. by 1.5%. In particular, if we also impose the strong 
constraint that fim/i^ is a constant, this model would 
have a value of H.^ that would be slightly lower than the 
true value, namely ~ 0.735. This would, in turn, make 
us overestimate the value of dA{z) to all redshifts (see 
Fig. 6). 

Now consider the BAO ratio at some redshift Zbao- 
For the sound horizon at the drag epoch, we would 



again make an overestimate by the same factor, ^ 1.5%. 
But at the same time, we overestimate dA{zBAo) or 
c/ H{zbao) by an almost similar factor because $7^ has 
shifted (see Figs. 6 & Figs. 7). Therefore, the trans- 
verse ratio rs{zD)/dA{zBAo) or the line-of-sight ratio 
rs{zD)H{zBAo) estimated in the wrong model will be 
rather close to the true values, thereby making it hard to 
detect the presence of early dark energy. In fact, these 
figures show that deviations from the fiducial ratio occur 
only at the 0.2% level at z < 0.5. Because the angular 
scale of the CMB acoustic peak is a precisely measured 
quantity, it is expected that even if the strong CMB prior 
is relaxed, a MCMC type exploration of the parameter 
space would reveal a similar insensitivity of the BAO to 

ni [16]. 

For the lensing ratio, the situation is even worse, as 
depicted in the right panel of Fig. 7. Here, due to over- 
all shift of all distance scales, the difference between the 
inferred and the true model is less than 0.1% at all lens 
redshifts. 
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FIG. 6: Left: Comoving angular diameter distance at the high redshift end for various models. The solid curve corresponds 
to the fiducial early dark energy model with fjjj — 0.742, wo = — 1 and SIJ^ — 0.03. The horizontal line indicates the value of 
dAizo) at the last scattering surface required by a wrong model with = 0, to keep the CMB angular scale, 9a, constant. The 
dotted line represents a model with all parameters kept same as the fiducial model, except which is set to zero. This falls 
short of the required dA{zo) and hence the only free parameter in the model, f2°, has to be adjusted to achieve the required 
dA{zo)- The final model that would be wrongly inferred by matching the CMB acoustic scale, has = 0.735 and is shown by 
the dashed line. Right: Fractional difference in the comoving angular diameter distance dA{z) from the fiducial model. The 
dotted line represents the fractional error for the same model as shown by the dotted line on the left plot. As expected the 
difference in this case goes to zero at low redshift. The dashed line shows fractional difference in the wrongly inferred model. 
Note that since this model had its 0° shifted low, it overestimates the true distances by ~ 1.5% for a > 1 and by 1.3% for 
z < 0.5. 
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FIG. 7: Left: Fractional difference from the fiducial of the Hubble scale c/H{z) in the wrongly inferred model of Fig. 6. Right: 
Same as left, for the lensing ratio as a function of the lens redshift. 



